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It is widely assumed that the observed reduction of the magnetic field of millisecond 
pulsars can be connected to the accretion phase during which the pulsar is spun up by 
mass accretion from a companion. A wide variety of reduction mechanisms have been 
proposed, including the burial of the field by a magnetic mountain, formed when the 
accreted matter is confined to the poles by the tension of the stellar magnetic field. A 
magnetic mountain effectively screens the magnetic dipole moment. On the other hand, 
observational data suggests that accreting neutron stars are sources of gravitational 
waves, and magnetic mountains are a natural source of a time-dependent quadrupole 
moment. We show that the emission is sufficiently strong to be detectable by current 
and next generation long-baseline interferometers. Preliminary results from fully three- 
dimensional magnetohydrodynamic (MHD) simulations are presented. We find that the 
initial axisymmetric state relaxes into a nearly axisymmetric configuration via toroidal 
OO , magnetic modes. A substantial quadrupole moment is still present in the final state, 

C ~ ~) ■ which is stable (in ideal MHD) yet highly distorted. 



1. Introduction 



Despite considerable effort, an unequivocal direct detection of gravitational waves 
(GW) is yet to be achieved. The expected wave strain is several orders of magnitude 
weaker than the sensitivity of current interferometric detectors. 1 One possibility is 
to coherently integrate the signal of a continuous source. In this case, the signal- 
to-noise ratio increases with the square root of the observation time. 2 A variety 
of physical mechanisms for the generation of continuous gravitational waves have 
been proposed, 3 ' 4 among them nonaxisymmetric distortions of the neutron star 
crust, either due to temperature variations 5 ' 6 or strong magnetic fields, r-mode 
instabilities, 8-11 or free precession. 12 ' 13 

A promising GW source was recently suggested by two of us. 14 Matter accreting 
onto a neutron star in a low-mass X-ray binary (LMXB) accumulates at the mag- 
netic poles until the latitudinal pressure gradient overcomes the magnetic tension 
and the plasma spreads equatorwards. The frozen-in magnetic field is carried along 
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with the spreading matter, and is therefore compressed, until the magnetic tension 
is again able to counterbalance the thermal pressure. This configuration is termed a 
magnetic mountain. 15 During this process, the magnetic dipole moment decreases 
with accreted mass, consistent with observational data. 15 ' 16 

We examine the prospect of detection of GW from magnetic mountains in section 
[2] In section[3j we present preliminary results from fully three-dimensional magneto- 
hydrodynamic (MHD) simulations to test the mountains stability. 

2. Gravitational radiation 

A typical mountain with M a 1CP 4 M0 and pre-accretion dipolar magnetic field 
of B = 10 12 G can provide 14 a gravitational ellipticity e = \h — h\/h ~ 10~ 5 , 
where I\ and I3 are the principal moments of inerta. e is considerably higher than 
the deformation a conventional neutron star could sustain (e k, 1CP 7 ) via its free 
elastic response and is only surpassed by exotic solid strange stars 17 (e ~ 10~ 4 ). 

The characteristic GW strain 18 is defined as h c = (128tt 4 / lb^^GI zz j ' 2 e/(dc 4 ), 
where I zz « 10 45 g cm 2 is the principal moment of inertia, / the spin frequency, 
and d the distance to the object. Fig. [1] (left) shows h c as a function of / at a 
distance of d = lkpc and a mountain mass of 1O _8 M0 < M a < 1O _2 M0 together 
with the design sensitivities of LIGO and advanced LIGO for a coherent integration 
time@ of 10 7 s. A mountain with M a « 1O~ 6 M should be detectable by LIGO for 
/ > 200 Hz. This is also consistent with the observed cutoff at ~ 700 Hz in the spin 
frequency distribution of LMXBs — much slower than the breakup frequency. 5 

3. Three-dimensional hydromagnetic stability 

Surprisingly, the distorted magnetic configuration is stable to axisymmetric 
modes. 19 However, the full three-dimensional stability is yet to be examined. We 
perform three-dimensional simulations by loading the axisymmetric configuration 
into the ideal MHD-code zeus-mp. A preliminary result is displayed in Fig. [T] (right). 
Shown is the time evolution of the three cartesian quadrupole moments Q22 , Q33 , 
and Q12, defined as Qij = f d 3 x' (3x' i x'j — r' 2 5ij)p(x'), where p denotes the den- 
sity. After readjusting initially, the system settles down into a state that still has a 
sizeable quadrupole moment. Furthermore, the small magnitude of the off-diagonal 
element Qi 2 suggests that the mountain is still nearly (within ml %) symmetric 
about the magnetic axis. 

We tentatively interpret these results as a preliminary proof of three-dimensional 
stability. However, the influence of resistivity still needs to be examined. Resistive 
ballooning and resistive Rayleigh- Taylor modes may allow plasma slippage on a 
short timescale. 20 Non-ideal MHD simulations to investigate these effects are cur- 
rently under way. 



a This is currently too optimistic due to computational limitations. The S2 run managed to reduce 
a five hour stream of data. 3 Improvements are expected using added computational resources and 
hierarchical search strategies. 
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(a) 

Fig. 1. (a) Characteristic wave strain h c for mountain masses 10- 8 M Q < M a < 10- 2 M Q along 
with the design sensitivity of LIGO and Advanced LIGO for a 99 per cent confidence and an 
integration time of 10 7 s. (b) Time evolution of the quadrupole tensor. The left and righ axes 
scale the diagonal and off-diagonal components, respectively. The time base is the radial Alfven 
crossing time t\ = 5.4 X 10 — 7 s. 
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